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IONOSPHERIC LIMITATIONS TO TIME TRANSFER BY SATELLITE 
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ABSTRACT 

The ionosphere can contribute appreciable group 
delay and pnase cnange to rauic signals traversing 
it; this can constitute a fundamental limitation 
to tne accuracy ot time and rrequency measurements 
using satellites. Because of the dispersive nature 
of the ionosphere, the amount ot delay is strongly 
frequency-dependent (l/f2j, the 1.5 GHz 

frequency band used in the GPS systems the 
vertical-incidence ionospheric excess delay is 
typically 5 nanoseconds during the daytime (based 
on a total electron content ot 10 ^/ m 2 ) > and 
20>fc of that at night. Even at X-band,the total 
daytime ionospheric delay of about 1 nanosecond 
is enough to make compensation necessary 'or 
implementation ot extremely precise time transfer 
schemes such as the coherent satellite link 
proposed by Knowles. Calibration using models is an 
unreliable procedure because of the variable nature of 
tne ionosphere. It is possible in principle to provide 
a self-calibration by observing at two frequencies 
simultaneously. While this tecnnique has on 
occasion proven successful in reducing ionospheric 
errors, a tuller understanding ot the underlying 
phenomenon is necessary in order to understand 
the basic limitations it mandates in time transfer 
techniques . 

The ionosphere is known to be highly variable, 
both in terms ot large-scale changes and in 
terms ot smaller-scale irregularities. Much 
work in recent years has focused on 
pathological 1 irregularities in the polar 
regions and equatorial region that amount to 
variations ol 1U)6 to several hundred percent 
in the total electron content. Recent 


229 




} 

t 

? 

i 


I 

i 


* 



efforts to investigate ionospheric disturbances 
using radio astronomy interf eromters have shown 
that the mid-latitude region, where most time 
transfer experiments take place, also has a 
prevalent irregularity distribution. This work 
has been undertaken by our NKL, group using 
interferometers at Green bank, w. Vh, and 
Magdalena, NM, and also by a Dutch group using 
the Westerbork interferometer (3). These 
experiments measure the differential electron 
content Detween two or more antenna locations, 
and relate more closely to the ionospheric effect 
seen by a time transfer experiment than most 
other techniques, easily measurable ionospheric 
irregularities nave been seen in raid-latitude 
regions. Our group has observed a typical 
differential electron content of about 10l5/m2 
over a baseline of 35 km. A day-night correlation 
with TlC values is observed; the night-time 
ionosphere is more irregular during periods of 
nign geomagnetic activity. by combining results 
from our work with that ol the Dutch group, who 
used a baseline that was considerably shorter, 
it is possible to obtain information about the 
size distribution ot the irregularities. The 
irregularity amplitude is approximately 
proportional to the baseline over the range 
studied. Only limited information is available 
aDout another parameter of interest, tae time 
scale distribution ot these irregularites . 

Irregularities have been observed on all time 
scale* from seconds to hours, as illustrated by 
data samples. Ine typical irregularity in the 
mid-latitude region is about 1% of the total 
electron content. While this is less than 
the extreme variations seen in polar regions, 
calibration of these anomalies must also take 
place to achieve full accuracy in time transfer 
and its associated use for position location. 

INTRODUCTION 

The ionosphere can contribute appreciable group delay to signals tiaversing 

it. This delay depends on both the frequency and the state ot the ionosphere. 

The equation tor the refractive index of tne ionosphere can be expressed as: 

M 1 "« * *2 n e B cose + 0(l/f 4 ). (1) 
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where n fc is the electron density 

b is the earth's magnetic field 

0 is the angle between the wave normal ana the magnetic field, 
and 

^1» ^2 are numerical constants iMathur at al . , 197U). 

It n e is replaced by a typical total electron content along a vertical 
column (TtC), an estimate can be obtained of the additional group delay due to 
the ionosphere, buch a curve is plotted in figure 1 for a typical daytime TtC 
of 2 X 10^/m^. Tne curve is for a vertical path and is multiplied by the 
secant ot the zenith distance of the line-oi-sigftt . The dominant term from 
equation 11) has a shape inversely proportional to frequency 3quared, ana a 
total magnitude of about 5 nanoseconds at the Global Positioning bystem U : Pb) 
frequency. The typical night-time delay is about 25/o oi this. Kiobucnar 
(1984) discusses the limitations in attempting to correct for the delay with a 
model due to tne variations in TtC. buch a model can be expected to be 
accurate to about + 25*. 

A method often mentioned tor correcting accurately for the ionospheric group 
delay is dual-frequency measurements, because of the dispersive nature of the 
ionospheric delay, measurements at two frequencies can be useu to determine 
what the delay constant is and correct tor it. In principle this method can 
be very accurate, because it measures the TtC along tne same line ot sight 
that is to be used for the observation. Although the dual-frequency method 
has been used in a number of experiments, and is to be implemented for GPb, 
its accuracy has not been critically examined. Possible inaccuracies in this 
method include error multiplication caused by a finite frequency uirrerence, 
and the effect of higher-order terms in equation (1). Figure 1 contains an 
estimate ot the residual comoonent if th^s correction is accurate to une 

residual component ot concern is the l/f component m equation 11). 1ms 
component gives rise to two possible phase and group velocities in a wave 
propagating through tne ionosphere* Tne airterence in velocities depenus on 
the angle ot the propagation direction to the earth's magnetic field. This 
effect gives rise to Faraday polarization rotation, but it also causes a 
difference in group delay between the two modes. Figure 1 shows the 
"worst-case" magnitude ot tms ettect. Wnile the net time delay change caused 
is zero, a propagation bifurcation results that causes a dual arrival time 
measurement. inis car. oe a serious limitation on lower-frequency timing 
measurements . 

When considering the limitations provided oy the ionosphere to time transfer 
measurements, it is important to consider the complete temporal and spatial 
spectrum ot ionospheric irregularities. Klobuchar (1964) discusses the 
large-scale (global) component of TLC variations. Recent experimental 
information has been accumulating about medium-scale (lOu m - iOoO tcm) 
irregularities in TtC that are important in contributing to time measurement 
inaccuracies. It has been deduced from in-situ satellite measurements that 
medium-scale irregularities of 1-iO X of TLC pxist commonly at ail times in 
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tne arctic, and during nighttime in equatorial regions (Szuszczewicz et al, 
1963) . Although it has been commonly assumed that medium-scale 
irregularities m the "normal" mid-latitude ionosphere are insignificant , 
recent research using radio astronomy interferometers have shown that this is 
not the case. 


In a radio interferometer, two or more large radio antennas are pointed at the 
same distant natural radio source. The outputs of these antennas are combined 
coherently, after correction for differential delay, to make a Michaeison 
interferometer. The output of this interferometer measures any source of 
differential phase delay. Since the ionospnere is the dominant dispersive 
effect, it is possible to measure the difference in total electron contert by 
combining measurements at two frequencies in an appropriate manner. To this 
end, the Naval Observatory's Green bank, West Virginia interferometer and the 
newer Very Large Array facility in New Mexico have been used to conduct 
preliminary investigations into the feasibility and usefulness of such a 
tool. The radio interferometer system at Green tank, whose primary mission is 
to measure geodetic parameters, observes contimouslya succession of radio 
sources at quasi-random points n the local sky. IKiepczynski et al., 1979). 
Simultaneous observations at 2695 MHz (11.1 cm. wavelength) and 8085 MHz (3.7 
cm. wavelength) provide a measurement ot the differential electron content 
between the ionospheric paths ot any of the possible antenna pairs. The 
continuous avaiiabiliity of data from tuis array makes its use attractive for 
ionospheric monitoring. At the time tne present data was recorded, the array 
consisted of 3 antennas ot 85-ioot uiaraeter located along a 20 kilometer 
track that is oriented at an azimuth of 62.05 degrees, and one remote station 
with a 45-toot diameter antenna located about ob kiu away. An additional 
remote station has since been added. The Very Large Array, described in 
detail by Thompson et al (1980), is an array of 21 antennas located near 
bocorro, NM, and designed with maximum flexibility for a variety of 
as trophysieal experiments. its 27 antennas allow significantly more complete 
spatial sampling of ionospheric phase dilierences than the Green bank array. 
Tne maximum baseline available is approximately the same as the Green bank 
array • 

The deferential phase path between two antennas due to the ionosphere is: 
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so that the differential phase aelay is proportional to the difference in 
total electron content along the two parallel paths. In the context of this 
article, D.L.C. will be used to refer to the experimentally determined 
difference in ionospheric path length as between the locations ot different 
antennas. It may be due to a number of different mechanisms; e.g., 
ionospheric gradients (Komesaroff, I960), irregularities, diurnal effects on 
TEC, etc. This quantity \f^n e dl -f^ n e dl) will nenceforth be called 
differential electron content = D.L.C. Converting from measurement units of 
delay to units of phase, a power of frequency is cancelled, so that: 
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where K - 2 * M 
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and A 4> i is the measured phase change due to the ionosphere. 

Another important parameter to cousider is the efiective area intercepted by 
each antenna primary beam as it transits the ionosphere, as well as the 
overlap between beams. The ditferential phase eftects on an interferometer 
from the ionosphere are described by a special case ot scintillation- 
scattering-diffraction theory; in particular one may expect thin-screen 
diffraction theory to serve as a basis. This type of theory is complex and 
has been discussed by a number of authors l Crane, 1977; Kufenach, 19751, 
although relatively few l Whale and Gardiner, 19661 discuss the expected pnase 
perturbations. While it is beyond the scope ot this article to develop a 
complete theory tor this phenomenon, the basic remark can be made tnat if the 
two antenna beams do not overlap the phase fluctuations will be fully 
developed, while for beams that largely overlap the phase fluctuations will be 
largely cancelled. With the 27 antennas of the VLA, it is possible to do a 
more complete reconstruction of the ionosphere. If the antenna oeams do not 
intersect while passing through the ionosphere, 27 independent points in the 
ionosphere, located along a Y-contigurat ion at distances from each other of up 
to 35 km, will be sampled by the VLA. 

The differential electron content measured by an interferometer pair is a sum 
of the true difference in the free electron density between the paths over 
each antenna, and a second order difference in between the paths over each 
antenna, and a second order difference in the path length caused by the 
curvature of the earth. The formula for the curvature effect may be shown to 
be in phase units 

A<*»= 4.85 x lO” 11 (TF.C) . / 1 - 1 \ (4) 

f lsin0 2 sinGj) 

where 0^ 02 are the elevation ot the source at each antenna^ t is the 
frequency in MHz, and A <t> is the differential phase in units of degrees. The 
difference tl/sin0 2 - i/sin 0^) can be either positive or negative, 
depending on the source geometry, but is always less than d/r where d is the 
distance between antennas, and r is the earth's radius, bince this is a 
refraction effect due to the sphericity of tne earth, it does not depend on 
the height of the electron layers doing the refracting. The total measured 
differential electron content is thus the sum ot the differential electron 
content due to variations in total electron content between the two locations 
and that due to the sphericity ot the earth. 

For the baselines used in our data sample, both effects are of comparable 
magnitude and must be considered. To measure the irregularity component, the 
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other must be subtracted. In order to do this, some estimate must be made ot 
the TEC. In analyzing the Green bank data TEC values obtained trom satellite 
Faraday rotation data by the Air Force Global Weather Control (1961) were used 
for this purpose. Values obtained at bagamore hill, HA were used as the 
closest approximation available to the Green bank, WV site. A time correction 
was made in the TEC data to compensate for the difference in longitude. In 
order to simplify the date analyses only phases trom two baselines, one of 33 
km and one of 1.3 km, were considered. Data taken with the U.b.N.O.'s 
Greenbank interferometer during November 1963 together witn smaller samples 
from 1980-81, was analyzed to isolate ionospheric effects, lhe periods 
include a wide range of geomagnetic activity. 

Figure 2 shows the estimates of variance of estimates of K for the long and 
short Green bank baseline tor the month ot November 1983. Data has been 
grouped into periods ol morning twilight, day, evening twilight, and night to 
emphasize diurnal effects. Occasional data points have been bridged. Both 
the higher absolute value of the variance on the long baseline and its 
increased variability are evident. Figure 3 shows a weighted average of the 
variance foi *:he four time periods tor the entire month. The day variance for 
the long baseline is clearly greater during the daytime, while this eftect is 
not present for the short baseline* An obvious interpretation ot this is that 
D.E.C. *s of 5xl0l^/m 2 are consistently present during the daytime over a 
35 km baseline, while irregularities over a 1.3 km baseline are consistently 
less than 1.5xl()14/ m 2 # xne phase or the erfect agrees with that of TEC 
measurements, which is normally four times higher during daylight hours. 
Although the constant portion of the long baseline variance may be due to the 
ionosphere, it is more reasonable to ascribe this to equipment effects. 

Althougn the variance correlates well with the diurnal TEC effect, it would be 
expected to correlate more directly with an index that would predict 
irregularities in the ionosphere. In Figure 4 the daytime long baseline 
variance is plotted on tne same graph as the geomagnetic index Ap. a 
general correlation is evident. No noticeable correlation of the daytime 
variance with the daytime maximum TEc was evident. This is an indication that 
irregularities are being measured, rather than some differential of the normal 
daily gradient. 

Another method of displaying tne results of the experiment is to follow the 
tracks of one or more sources and compare the measured differential 
ionospheric eftect witn that predicted trom independent total electron content 
data. This is done in Figures 5 and 6 for the radio source 0355+5U8, which 
was observed at night, and tne radio source 1749+701, whicn was observed 
during the day. Also plotted on these figures is the differential path length 
ionospheric effect predicted from the total electron content data. The data 
from the nighttime source is seen to be generally smoother, although variable 
*rora night to night. That from the daytime source always, and that from the 
nighttime source sometimes, has fluctuations of 1-2x10^5 £ n differential 
electron content. The data is too sparsely sampled to enable a detailed 
Fourier analysis of the time scale of the fluctuations, but it appears that a 
typical time scale is of the order of 1-2 hours. For the baseline of 35 km, 
the observed 1-2 hour time scale means that tne ionospheric irregularities 
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observed are typically either changing with that time scale or moving with a 
velocity of at least 20 km/hr (7m/sec) with respect to the earth's surface. 
This estimated speed is much lower than currently accepted values tor the 
ionospheric superrotation rate ot 50-150 m/sec. 

A sample detection of an ionospheric irregularity using the VLA is shown in 
Figure 7. A method similar to tnat described above has been used to separate 
phase excursions at frequencies of 1.4 GHz and 5 GHz into dispersive, or 
ionospheric, and non-dispersive, or tropospheric terms. The uata shown 
includes spacings of 17 km., 9 km. , 3 km., and 0.4 km., during a period of 24 
hours beginning at 08 hours local time on 21 January 1979. The figure snows 
observed phase in units of degrees at 1.4 uhz (.100 degrees = 2.5x10^^ 

D.L.C.). During this period, the source 1311+678, which is circumpolar, was 
observed until about 14.7 hours local time; then source 0352+398 was observed 
until the first source was again visible at about 0.7 hours. Source-change 
times are indicated by the long, vertical bars extending through tne figure. 

The deep protrusion in the 17 km. and 9 km. ionospheric outputs may be 
interpreted as a gradient in TtC which lasted for a duration oi about 7 hours, 
with a mean epoch of about 16 hours local time. Its apparent intensity is 
magnified by the fact that the observations at that time were low in the sky. 
After correcting for a secant effect multiplication ot about 3, an 
irregularity of 6x11)14 i n u.E. c. is indicated; this is consistent with the 
Green bank results. 

The 9 km spacing is generally similar to the 17 km spacing, but is decreased 
in amplitude by about 30)4, as would be expected for an ionospheric gradient. 
Thus, the total D.L.G. may be larger depending on the size ot the 
irregularitv. A consistency check is provided by the tropospheric output, 
which snows no significant change during this period. 

OISCUbSION 

fhe radio astronomy technique, can measure ionospheric structure on a variety 
of scales, and provides a measure of an effect integrated through the 
ionosphere, rather than the marginal electron content measurement made by 
either bottomside or topside (keinisch and Xuegin, 1982) sounders. During the 
observational periods described here, irregularities ot the order of from 
2x10^ to 2.5x1015 in differential electron content were observed. This 
amounts to about 0.5 to 3 i ot the measured ThC, and indicates the frequent 
presence ot significant irregular structure in the ionosphere on this scale. 

A Dutch group (bpoeistia and Kelder, 1984) have recently used the Westerbork 
interferometer to measure ionospheric disturbances. Although their longest 
baseline available was only 2.7 km., they were able to easily measure 
ionospheric irregularities. They found the typical size to be consistently 
greater than their maximum baseline, which is m agreement with the present 
size estimates. Their conclusions about the greater prevalence of 
disturbances during the daytime are in general accord with ours, although the 
Dutch group did not provide numerical estimates ot differential electron 
content. 
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Bo In the Green bank results and the sample of data from the VLa show clear 
evidence of horizontal ionospheric gradients in TEC persisting tor periods of 
several hours, while the Dutch group ooserved periods of less than fitty 
minutes. The mid-latitude radio interferometer observations clearly are not 
reiated to the ionospheric wind, which is known to have velocities of 50-130 
kilometers per second (Regan and Rodriquez, 1981; kino and Livingston, 1982). 
Iney are most easily accounted for by tne phenomenon of traveling ionospheric 
disturbances, or HD's. This class of disturbances, as pointed out by Evans 
et al 11983), includes fluctuations with a wiue variety of periods, which, in 
spite of many observations with various techniques, are not well-defined. 

From the limited observations undertaken by our group and the Dutch . group, it 
is possible to make a rough estimate of the approximate amplitude vs. distance 
relationship ot these medium-scale irregularities. This is showu in figure 
8. A similar graph that gives typical time scales is not yet possible due to 
inadequate data. 

SUMMARY 

Ionospheric compensation is necessary for the most precise time transfer and 
frequency measurements, with a group delay accuracy better than 10 
nanoseconds. A priori modeling is not accurate to better than 23&. The 
dual-frequency compensation method nolds promise, but has not been rigorously 
experimentally tested. Irregularities in the ionosphere must be included m 
the compensation process. 
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Fig* 2. Scatter in ionospheric phase measurments for Green bank ; 

interferometer for the month of November, 1983. Dotted line marks \ 
short baseline; solid line marks long baseline. Data is grouped into , 
morning twilight, day, evening twilight and night. ; 
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Fig. 3 All data of Figure 2 grouped to show behavior of variance as a 
function of time of day. 
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Fig, 4 Comparison oi Green bank long baseline variance tor daytime with 
geomagnetic inuex Ap, ixorinalized units used on ordinate. 
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Fig. 5 Differential electron content measurement as a function of time for 
035!>+3U8, a nighttime source. Observations taken during early 
December, 1981. 
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Differential electron content measurement as a function of time tor 
1749+701, a daytime source. Period ot observations same as Fig. 5. 
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Fig* 7 Sample ionospheric measurements made using a series of baselines with 
the Very Large Array , illustrating the detection of an ionospheric 
disturbance. Tne graph shows both non-dispersive , or tropospheric 
and dispersive or ionospheric differential phases* Long vertical 
bars mark source changes* 
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Fig. 8 Amplitude vs. size parameters of ionospheric irregularities. 
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QUESTIONS AND ANSWERS 


LARRY D'ADDARIO, NATIONAL RADIO ASTRONOMY OBSERVATORY: This 
informati n about the variations between stations on relatively 
short baselines would also suggest that, for a single station, 
there would be a lot of variation in direction if you are looking 
at the sky. Thu3, knowledge of the zenith total electron content 
would not be a good predictor of the path delay in other 
directions where you are usually looking. Is that right? 

MR. KNOWLES: I would agree with that statement. 

MR. KLOBUCHAR: If I could interject a comment here, that's 
certainly true, but you have to look at what detail is true. If 
Steve would show the 1 st view graph.... (the view graph was 
shown here). You will notice that the two real points that you 
said were actually observed, were down around ten to the 
fifteenth. The point way up at the top here at around a thousand 
kilometers above ten to the seventeenth is actually a difference 
in the diurnal curve between two points. The total time delay 
might be on the order of ten to the seventeenth, and you are 
seeing effects here on the order of ten to the fifteenth, one or 
two percent. Ten to the sixteenth is half a nanosecond at L-band, 
so we are talking about much less than a nanosecond, differen- 
tial. If that's a problem for your time transfer, then you are 
in trouble, because you always get fluctuations of this order. 
For the VLBI folks, it is a serious problem, and Steve's work is 
very important in this regard, because this size of irregularity 
is one that really hasn't been measured before. It is very 
important to continue that kind of work, because there are 
irregularities at all levels, and no one up until this work has 
done anything in that distance range that I am aware of. The 
Dutch work and this work are the only ones. 
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